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Abstract 

We examine the possible extension of the parameter space of the minimal supersymmetric ex- 
tension of the Standard Model (MSSM), as expressed via the renormalization-group equations in 
terms of universal soft supersymmetry-breaking terms at the unification scale, to include tachyonic 
input scalar masses. Many models with negative masses-squared for scalars at the unification scale 
may be viable because the small sizes of the masses-squared allow them to change signs during 
the renormalization-group evolution to the electroweak scale. However, in many cases, there is, 
in addition to the electroweak vacuum, a much deeper high-scale vacuum located along some F- 
and D-fiat direction in the effective potential for the MSSM. We perform a numerical search for 
such vacua in both the CMSSM and the NUHM. We discuss the circumstances under which the 
existence of such a deep charge- and color-breaking vacuum is consistent with standard cosmology. 
A crucial role is played by the infiation-induced scalar masses, whereas thermal effects are often 
irrelevant. 
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I. MOTIVATION 



Understanding the allowed parameter space in versions of the Minimal Supersymmet- 
ric extension of the Standard Model (MSSM) with Grand Unified Theory (GUT) inspired 
boundary conditions is a research programme that gains motivation at the onset of the LHC. 
In the constrained MSSM (CMSSM), in which gaugino and scalar masses are unified at the 
GUT scale, and its generalizations with non- universal boundary conditions, e.g., for the 
Higgs scalar masses (NUHM), the identity and relic density of the lightest supersymmetric 
particle (LSP) place strong constraints on the parameter space [])]. For example, in models 
with small values of the unified scalar mass, mo, and large values of the gaugino mass, mi/2, 
the LSP is typically the lighter spartner of the r lepton. This region of the MSSM param- 
eter space would predict charged dark matter, and hence would be excluded if i?-parity is 
conserved and the fi is stable. However, this corner of parameter space might be allowed 
if the gravitino is lighter than the fi, and becomes the LSP j^, as may occur in mSUGRA 
models [s^. In this case, the fi decays to the gravitino LSP, and is subject to important 
astrophysical constraints that do not exclude this region. If the gravitino is the LSP, the 
CMSSM parameter space may even extend to negative values of mg [4.]^. In this paper, we 
argue against excluding all of this region because of unreasoning tachyophobia. 

One can define an effective MSSM model by specifying its mass parameters at the weak 
scale. In models with relatively light squark masses, the renormalization-group equations 
(RGBs) for the scalar masses may then lead to m^iQo) = 0, for some value of Qo in the range 
M\Y < Qo < Mgut, with rn?{Q) < for Q > Qq- This raises the question whether there is a 
dangerous charge- and color-breaking (CCB) vacuum ^. The answer to this question depends 
on two factors: whether there are potentially large logarithmic corrections to the potential 
which are not absorbed in the running of mass parameters, and whether there are significant 
non-renormalizable terms in the effective potential. In general these vacua, determined by 
minimization of the tree-level scalar potential, occur with vacuum expectation values (vevs) 
of order ~ \m?\ <^ Q^. However, m^{Q) < and the existence of F- and D-flat directions 
in the MSSM leads to runaway to v{Q) ~ Q for Q > Qq, where logarithmic corrections are 
small. Such a CCB vacuum would exist if there were no non-renormalizable terms in the 

^ Negative scalar masses-squared also appear at the high-eneroi scale Q in a version of mirage mediation Q . 
^ For discussions concerning CCB vacua in the CMSSM, see [71]. 
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effective potential, but such terms are in general present, and their magnitudes determine 
where the runaway vev is stabilized and hence whether the existence and location of such a 
CCB vacuum can be calculated reliably. Thus, as argued in 81] , for certain parameter choices, 
vevs will indeed be generated with v{Q) ~ Q for Q > Qq, ensuring the existence of the CCB 
vacuum unless new physics between the GUT and electroweak scale is introduced. Here we 
will enlarge upon this idea and identify the flat directions that are primarily responsible in 
the CMSSM and NUHM. 

However, it is possible that such CCBs could be tolerated j^, if the Universe would have 
fallen naturally into our false electroweak (EW) vacuum as the cosmological temperature 
decreased, and if the lifetime of this vacuum for tunnelling into the true CCB vacuum is much 
longer than the present age of the Universe. Whether the EW vacuum is in fact preferred by 
cosmology depends, in particular, on the scalar masses-squared generated during inflation. 
If these masses-squared are positive and of the order of the square of the Hubble parameter, 
the 'more symmetric' EW vacuum is favored. On the other hand, if these are negative, as 
in the Affieck-Dine scenario for baryogenesis 10| , the Universe would remain trapped in the 
true CCB vacuum [l^. 

The present article delineates regions of parameter space for which high-scale CCB vacua 
are present. We study both the CMSSM and less constrained models with non-universal 
Higgs masses (NUHM) Q, Q. For each choice of GUT-scale parameters, we follow the 
analysis of l4] for determining the set of problematic fiat directions and the lowest order of 
GUT- or Planck-scale non-renormalizable operators whose appearance might lift the fiat di- 
rection. The existence of CCB vacua will depend on both the order of the non-renormalizable 
operator and the fundamental scale associated with it. In the case of the CMSSM, there 
are wedges of parameter space with rriQ < where no MSSM sfermion is tachyonic at the 
EW scale. However, these regions generally have calculable CCB vacua, assuming that the 
mass scale M^, in the non-renormalizable interaction that stabilizes the high-scale vacuum is 
greater than or equal to Mqut- Such regions would be acceptable in suitable cosmological 
scenarios that populate exclusively the EW vacuum. Similar questionable tachyonic regions 
occur also in the NUHM, even if ttiq > 0, due to the (independent) squared masses for 
Higgses being negative. However, we emphasize again that these regions would be accept- 
able in cosmological scenarios that avoid populating high-scale vacua: one should not be 
unreasoningly tachyophobic. 
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II. CCBS ALONG FLAT DIRECTIONS IN THE MSSM 



The tree-level scalar potential can be written schematically (including soft 
supersymmetry-breaking contributions) as 

V = \rh\Q)ct>' + \\{Q)ct>\ (2.1) 

where, for simplicity, we have neglected cubic terms, and A may vanish along some di- 
rections in (j) space. We assume as well that some set of soft supersymmetry-breaking 
scalar masses have fh?{Q) < for Q above some scale Qo; for a particular choice of GUT- 
scale parameters in the CMSSM or NUHM. Along a generic direction in field space, with 
\{Q) ~ g'^ where g is some gauge coupling, minimization of the tree- level potential yields 
a vev at = v{Q) = [—m'^{Q)/\{Q)Y^'^ which is of similar order to the supersymmetry- 
breaking scale, m3/2- However, one- loop corrections to the potential will have the form 
AV ~ rh"^ (Q) V (Q)^ \ii[v {QY /Q^] and, for v ^ Qo, the large logarithms may reduce signifi- 
cantly the reliability of the CCB vacuum calculation. 

On the other hand, the appearance of negative mass-squared would have important 
implications for the moduli rj of the flat directions in the MSSM where A = 0, since the 
tree-level solution yields a runaway vev. However, we expect high-scale non-renomalizable 
operators to regulate the runaway behaviour in such a case. In the presence of a non- 
renormalizable superpotential VFnr of degree n (see Appendix A), the effective potential 
becomes 

= + (2-2) 

and the runaway direction for the modulus f] is stabilized. In (12. 2p . u is an 0{1) dimensionless 
coupling constant and M^, is the cutoff scale associated with the dynamics that generates 
Wnr- Depending on the circumstance, it could be between the unification scale and the 
Planck scale. The modulus r] then acquires a vev whose order of magnitude is given by: 

V = [-m2^2{n-3)j 1^1) _ (2.3) 

General properties of fiat directions and the operators which lift them are found in Ap- 
pendix A, and examples of specific MSSM flat directions are given in Appendix B. 
The one-loop correction to the scalar potential can be written as 

A^i_ioop = -^g^v^ml^ In (2.4) 
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where is the number of multiplets that get masses at the scale m = gv. The cancellation 
between boson and fermion loops has been taken into account, as is apparent from the factor 
m^yg. This is an average soft supersymmetry-breaking mass-squared, determined from the 
one-loop contributions of the states. 

If V is small, the overall loop correction can be large relative to the tree-level poten- 
tial (12.21) . which satisfies V ~ (9(m^f^) near the minimum, rendering the tree-level analysis 
unreliable. We introduce a parameter e that represents the boundary where this large- 
logarithm problem occurs: 

V ^ eQ ^ unreliable, v ^ eQ ^ reliable. (2.5) 

Since v shuts off at low Q in the models we consider, it is possible in principle that the 
tree-level analysis may never be reliable. We have nothing to say about such models in this 
article, except that they are not excluded by our analysis. On the other hand, when we find 
a vev V and there is a regime of Q for which this tree-level result is reliable, a vacuum state 
really exists at some large value of the modulus 77. The existence of this deep vacuum at 
large field values — generally color or charge breaking (CCB) — has cosmological implications, 
which we discuss below. To determine the parameter e, we note that the loop correction 
(12. 4p is comparable to the tree term m^f^ when: 

'V/ln^~-l ^ e = exp(-327rVA^(7'). (2.6) 



When V > eQ, we can reliably state that at least two vacua exist: the electroweak (good) 
vacuum and the high-scale (bad) one. 

Now a word on numbers: in (12.61) . we should take 1/10 < < 1/2, corresponding to 
the running constants between the EW and GUT scales. The number N depends on the 
flat direction, but would typically range between (9(10) and (9(100). For the larger value of 
A^, we have e ^ O{10~^); moreover, e falls exponentially as decreases, so it is generically 
much smaller than 10"'^. In other words, the loop factor suppression g'^/6A7c'^ means that 
the logarithmic enhancement must be quite large before the tree-level analysis becomes 
unreliable. As a consequence, the RG-improved tree-level analysis is generally a reliable 
indicator of the high-scale vacuum. 

As an example, consider the LLe'^ flat directions in the MSSM. Eleven chiral multiplets 
participate in the mass matrix once renormalizable superpotential terms are accounted for, 



namely e^, z/,;, e?, H^, H7 . The moduli space is once D- and F-flatness constraints are 
taken into account J^]. Thus eight chiral multiplets get 0{v) masses, although in some 
cases they are suppressed by very small Yukawa couplings y'^, corresponding to D-moduli 
whose flatness is lifted by the renormalizable superpotential ^. Furthermore, SU{2) x f/(l) 
gauge multiplets get 0{gv) masses. Each chiral multiplet contributes 2 to A^, because it 
includes a complex scalar and a Weyl fermion. A similar contribution comes from each of 
the four vector multiplets contained in SU(2) x U{1). Thus = 2(8 + 4) = 24. Even if we 
assume the GUT-scale value = 1/2, we already get e ~ 4 x 10^^^. In actuality, for many 
of the masses we should use (y^)^ <^ and e is even much smaller still. The result is that 
we trust the tree-level analysis for practically all v/Q: it is a robust result that we have at 
least these two vacua, and the cosmological arguments apply. 

Fig. [1] represents schematically several generic cases for the possible variation with Q of 
the vev v{Q). For each of the curves #1 - #4, v{Q) 7^ at large Q, but v{Q) ^ at small 
Q. In each of the cases #1 - #3, the condition (12.51) is satisfied over some range of Q and 
hence, according to the arguments presented above, the existence of a high-scale vacuum 
can be predicted reliably in each of these cases. These are representative of what happens in 
much of the parameter space that we study below in specific CMSSM and NUHM scenarios. 
On the other hand, when we consider curve #4 in Fig. [1], we see that, over the entire range of 
Q, the tree-level prediction cannot be trusted, owing to the persistence of large logarithms. 
For this reason, we exercise caution and choose not to apply the cosmological constraint 
in such a case. In some specific cases, a one- or two-loop analysis might be reliable, and 
the possibility of a high-scale (bad) minimum could be examined in more detail. However, 
studies of such possibilities lie beyond the scope of this work, where we restrict ourselves 
conservatively to the criterion (12. 5p . 



III. THE COSMOLOGICAL CONSTRAINT 



Having established the criteria for determining the existence of a high-scale vacuum 
(which we apply to specific CMSSM and NUHM models in the following two sections), 
we now turn our attention to the cosmology of models with high-scale vacua, asking how 

The LLe moduli that inhabit C"^ arc Hftcd by the non-renormahzable superpotential and get 0{rh) masses, 
which can be ignored in the loop corrections. 
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FIG. 1: Curves of the flat direction vev v = v{Q) (solid) versus various constraint curves (dashed). 
In regions where v{Q) lies above the dashed lines, the tree-level effective potential is a reliable 
indicator for the existence of a high-scale (bad) vacuum, according to the criterion \2. 5|) . 



problematic their existence may be. Such a global CCB vacuum with f ~ Q ^ My/ is 
separated from the local charge- and color- conserving minimum at the origin (or the EW 
scale if the Higgses are the scalars in question) by a potential barrier. As argued in 9|, this 
model remains perfectly acceptable if the Universe is trapped in the local minimum near or 
at the origin, as the timescale for producing a bubble of lower high-scale vacuum is generally 
much longer than the age of the Universe. 

However, in models with negative mo(MGUT), there are many such vacua corresponding 
to different flat directions, and each has a larger domain of attraction than that of the 
EW vacuum. Thus, for arbitrary initial conditions, the system is overwhelmingly likely 
to end up in one of the 'bad' vacua. However, this question must be analyzed in some 
appropriate cosmological setting, which introduces two extra ingredients affecting the shape 
of the effective potential: inflation and thermal effects. During inflation, supersymmetry 
is broken by the vacuum energy, which results in an extra contribution to the soft scalar 
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masses, of the general form 

Am^ = cH^ , (3.1) 

where H is the Hubble expansion rate during the inflationary epoch. The dynamics of the 
flat directions depend crucially on the sign of c, which is model-dependent. For a minimal 
Kahler potential, c = 3. On the other hand, for a no-scale Kahler potential, the induced 
scalar masses are zero at the tree level and loop corrections generate c < for fiat directions 
not involving the stop fl?! . 

We now consider different possibilities for the coefficient c, considering first the possibility 
c > 0. We assume also that the initial value of the flat-direction field (p is 10^'' GeV, in which 
case V{(f)) -C Mnflaton, and the evolution of </> does not affect the Hubble constant significantly 
until the late stage when H ~ "^3/2. In this case, obeys the equation of motion 

4> + 3H(}) + V'{4>) = , (3.2) 

where if is a slowly-varying function of time that can be described by an adiabatic approx- 
imation. For our purposes, V' can be approximated by cH'^cf) until H becomes comparable 
to m^i2- The effect of the non-renormalizable term in V{(f)) is less significant. The general 
solution of (13. 2p is 



(pit) = Ci exp 



\h{-3 + 79^) t 



+ C2 exp 



\h{-3 - V9 - 4c) t 



(3.3) 



where C\^2 are determined by the initial conditions. 

We consider first the classical evolution. For c = 0(1), the magnitude of (p scales as 
(p{t) ~ exp[— (9(l)if)f:]. Thus, within 5-10 Hubble times, (p will be of order H even if its initial 
value was very large ^ and, within the next ~30 Hubble times, (p will be of the electroweak 
scale. 

We next take into account the de Sitter quantum fluctuations, which play an important 



role in the dynamics of flat directions [11]. For a scalar field of mass ~ if 15 1 



~ H"^. (3.4) 

It is then clear that at large 0, the classical damping dominates and when cp becomes 
comparable to if, these quantum fluctuations are as important as the classical evolution. 



Here we take H ~ 10"-i4 Qgy ^^^^ ^ i^n q^^. 
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From that moment on, undergoes random oscillations of order H per Hubble time plus 
classical damping which decreases its magnitude by a factor 0{1) per Hubble time. As H{t) 
decreases, so does the amplitude of random oscillations about the origin. At the time when 
the soft term rriQ becomes relevant (which is after the 'infiaton oscillation' era), the amplitude 
of oscillations is O {171^/2) ■ The field settles at the origin and the quantum oscillations are too 
small to reach the barrier separating the two minima, which is further than O {1713/2) away 
from the origin. Therefore, the presence of a deep CCB minimum would not be problematic 
for c ~ 1. 

However, these conclusions do not in general apply for small positive c: < c <^ 1, which 
is a borderline case. First, the classical evolution is slow and, secondly, the amplitude of 
quantum oscillations is larger: (0^) ^ SH"^ /{Stt'^c). Whether the field settles at the origin 
at the end of inflation depends on further specifics of the inflationary model as well as the 
magnitude of c. 

For c < (see [l6^), the minimum of the potential during (as well as after) inflation 
is at large (p. Classical evolution will drive (p towards this minimum, whose position is a 
slowly-evolving function of time. At ~ ^^3/2, the field freezes at the CCB minimum, and 
the quantum fluctuations do not play any significant role. It is important to note also that 
thermal effects are irrelevant at large and cannot destabilize the CCB minimum. This 
is because all the fields couple and receive masses of order (multiplied by the gauge 
or Yukawa couplings). Thus they are heavy and cannot be thermalized. Consequently, no 
thermal mass term T^0^ is generated by thermal loops ^. 

We thus conclude that, for c < 0, the presence of deep CCB minima is ruled out by 
cosmological considerations. This applies in particular to models with the Heisenberg sym- 



metry 



18| . including no-scale models of supergravity 19|. Furthermore, this precludes the 



possibility of the Affleck-Dine (AD) mechanism for baryogenesis [10|, which requires nega- 



tive c [16|. Concretely, this applies to the cases #1 - #3 of Fig. [H in which the existence of 
such a bad high-scale CCB can be predicted reliably. On the other hand, even these cases 
would not be excluded for positive c ~ 1, and possibly also for small positive c. 

We remark finally that we have neglected the trilinear soft supersymmetry-breaking A- 
terms in the above considerations. In general, one expects A-terms of order H to be gener- 



Note that, although (p itself is light, it has no self-interactions at the renormalizable level. 
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ated during inflation. If their magnitude is a few times larger than H, an additional local 
minimum at large appears even for positive c ~ 1. Given a large initial value of 0, the 
fleld will evolve to this minimum and remain there after inflation. Thus, deep CCB vacua 
can be problematic even for c ~ 1, if the A-terms are sufficiently large. 

In what follows, we will examine specific CMSSM and NUHM parameter sets for which 
< at some high renormalization scale, with the aim of elucidating whether one should 
worry about them. 

IV. TACHYONS IN THE CMSSM 

We begin by considering the parameter space of the CMSSM, in which all gaugino masses 
are unified at a common scale (where gauge-coupling unification occurs) with the common 
value mi/2- Similarly, all soft scalar masses are unified at the same renormalization scale 
with a common value mo, as are the trilinear terms with a common value Aq. The remaining 
free parameters are the ratio of the two Higgs vevs, parametrized by tan (3, and the sign of 
the fi parameter. 

For each choice of these four input parameters (plus the sign of fi), the low-energy spec- 
trum can be determined and compared with phenomenological and cosmological constraints. 
In Fig. [21 we show the {mi/2,mo) plane for two fixed values of tan (3, both for /i > and 
Aq = 0. The sign displayed on the vertical axis is actually the sign of ml, so what is displayed 
is ml/\mo\, strictly speaking. In panel (a), we have fixed tan/? = 10 and, for mg > 0, we see 
results common in many CMSSM studies |lj]. The dark (brown) shaded region corresponds 
to the region for which the stau is the LSP and as such this region is normally excluded 
unless the gravitino is in fact the LSP ^. The medium (green) shaded region, at low mi/2 is 
excluded by the constraint arising from the branching ratio of 6 — > 57. The vertical dashed 
line is the chargino mass contour at 104 GeV, and the nearly vertical dot-dashed line is the 
Higgs mass contour at 114 GeV, as obtained using FeynHiggs [21j]. Only regions to the right 
of these lines are allowed by LEP. The pink band bordered by black solid curves is the region 
where supersymmetric corrections to the Standard- Mo del calculation of {g — 2)^ match the 
experimental measurement of {g — 2)^ within 2-a uncertainties (between the dashed curves 

^ Note that, even in this case, much of the region shown is excluded due to effects during and after Big-Bang 
nucleosynthesis involving the bound state of the stau and He: see and references therein. 
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FIG. 2: Portions of the CMSSM {mi/2,TnQ) planes for (a) tan/? = 10 and (h) tan/3 = 50. The 
notations are described in the text. The lower pink region is excluded because the lighter stau would 
he tachyonic at the EW scale. The regions below the dash-dotted curves have problematic high-scale 
tachyons for the indicated values of if e = 1, while for e ^ 10^^ all regions with rrig < are 
excluded. 



agreement occurs at the l-cr level). Finally, in the (turquoise) shaded region that tracks the 
stau LSP boundary at large mi/2, the relic density of the lightest neutralino would lie in the 



221, if 



range of the cold dark matter density determined by WMAP and other observations 
this neutralino were the LSP. 

Also shown in Fig. [2] (a) is a large (pink) shaded region at low mi/2 and negative mg where 
one of the sfermions is tachyonic at the electroweak scale. This region is also excluded. Of 
particular interest to us here is the region where ttiq < but the shading indicates that the 
lighter stau is the lightest spartner of a Standard Model particle but is not tachyonic at 
the electroweak scale: we repeat that parts of this region may in principle be viable ^ if 
the gravitino is in fact the LSP. We perform our scan of the parameter space in this region 
looking for flat directions which could lead to a bad high-scale minimum of the potential. 

According to the analysis of [l4], the QQuue flat direction is not fully lifted until degree 
n = 9, so from Eq. (12.31) we see that v along this direction is quite close to the non- 
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renormalizable mass scale M^,. The dot-dashed curves in Fig. [2] labeled by M^, demarcate 
regions for which we find solutions to v{Q) = Q: for a given value of M^,, all regions below 
the corresponding curve admit solutions to v{Q) = Q. Thus, below the curves we find a bad 
CCB vacuum, that might invalidate the corresponding choice of parameters, subject to the 
cosmological considerations discussed in the previous Section. When M^, = Mp, all regions 
with TTiQ < have this problem. Because all of the scalars have mf < at high enough 
scale when ttIq < 0, this fiat direction always has a reliable bad CCB minimum when we 
take e < 10~^ for all > Mqut- Thus all CMSSM models with < and e < 10"^ are 
problematic according to our analysis. 

In panel (b) of Figure [2], we show the analogous parameter plane when tan/? = 50, /i > 
and Aq = 0. The part of the (mi/2, ""^o) plane shown is dominated by regions where the stau 
is the lightest sparticle apart from the gravitino (shaded dark brown) and where the stau is 
tachyonic at the EW scale (shaded pink). The standard CMSSM phenomenologically and 
cosmologically acceptable regions occur at mi/2 > 200 GeV, and so are not visible in the 
part of the (^1/2,^^0) plane displayed. The wedge-shaped brown region with mg < has a 
calculable high-scale CCB minimum, as was the case tan /? = 10. For this reason, we do not 
display planes with intermediate choices of tan jS. We have also scanned the CMSSM planes 
with Aq 7^ 0, and the results are qualitatively similar to those shown here. 

To summarize, because in practice e ^ 10"'^ for any possible fiat direction, and more 
especially for the QQuue one, the CMSSM with ml < always has a bad CCB vacuum. 
However, this may not be populated, when cosmological considerations are taken into ac- 
count. 

V. TACHYONS IN THE NUHM 

One alternative to the CMSSM is the NUHM, where the scalar partners of the quarks and 
leptons still unify at Mquti but the soft supersymmetry-breaking scalar masses associated 
with the two Higgs doublets do not. This class of models has, effectively, two additional free 
parameters relative to the CMSSM. These are often chosen to be the weak-scale values of /i 
and the Higgs pseudoscalar mass, mA- Whilst it is certainly possible within the context of 
the NUHM to choose < (leading to the same difficulty with the QQuue fiat direction 
as in the CMSSM), problems with CCB can already occur for certain choices of /i and itla 
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even when mg > 0. This is because, when the weak scale \fi\ is large, typically one or both 
of the Higgs squared masses ml 2 = /^^ + 2 is negative at the GUT scale, a problem that 
is accentuated at small rriA- In this case, the squark and slepton masses- squared remain 
positive throughout the RGE evolution, avoiding the QQuue flat direction problem, and 
vevs V <^ Q for all Q may develop along the H1H2 flat direction, because this is lifted by 
a lower-order non-renormalizable term (n = 4). As a consequence, a larger fraction of the 
parameter space is allowed. 

At scales Q not too far above the EW scale, one begins to see the EW vacuum. The 
negative masses-squared of MSSM Higgses can lead to unacceptable vacua with CCB also 
at this low scale. We do not include the details of such low-scale CCB vacua in the analysis, 
as this would require a much more careful treatment of loop corrections and contributions 
from all soft supersymmetry-breaking operators. For this reason, we terminate the NUHM 
analysis at Q = 10 TeV in our numerical studies. 

In Fig. [HI we show examples of {fi,mA) planes in the NUHM for fixed mi/2 = 300 GeV, 
mo = 100 GeV, and Aq = with tan/3 = 10 (panels a, c and e) and tan/5 = 50 (panels b, 
d and f). We use the same shadings as used for the CMSSM to denote regions excluded by 
b ^ s'-f (shaded medium (green)), which excludes much of the parameter space at yU < 0, 
and a region in which the lighter stau would be the lightest spartner of a Standard Model 
particle (shaded dark (brown)). This includes two areas with relatively small fi and rriA, 
for tan/5 = 10, and most of the left side of the plane, for tan/? = 50. New to this figure 
are regions shaded dark (blue) for which the sneutrino is the lightest spartner of a Standard 
Model particle, that are seen at large fi and rriA, for tan/3 = 10, and in the upper right 
of the planes, for tan/? = 50. At least parts of the regions with a light stau or sneutrino 
could be allowed if the gravitino is the LSP. Once again, the dashed vertical lines at small 
fi show the 104 GeV chargino mass contour. The thin blue lines show the contour where 
rriA = 2m^ and the regions of good neutralino relic density near these lines correspond to the 
rapid-annihilation funnel region familiar from the CMSSM at large tan (3. Other strips with 
an acceptable neutralino LSP relic density appear in the stau and sneutrino coannihilation 
regions, running parallel to boundaries of the brown and blue regions, and in a 'crossover' 
strip close to the chargino exclusion, where the neutralino is a mixed gaugino/higgsino state. 

We show in each panel of Fig. [3] three dot-dashed contours with differing M^,. The inner 
curves (with the lowest values of |/i|) correspond to = Mp, whereas the middle and 
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tan p = 10 , mi,2 = 300 , mp = 100, e = 1 tan (3 = 50 , mi/2 = 300 , mp = 100, e = 1 




II (GeV) 1^ (GeV) 



tan p = 10 , mi/2 = 300 , mp = 100, e = lO"* tan |3 = 50 , = 300 , mp = 100, e = lO"* 




Pi (GeV) \i (GeV) 



FIG. 3: Portions of the NUHM {fi,mA) planes for for mi/2 = 300 GeV, mo = 100 GeV, Aq = 

and (a, c, e) tan/? = 10, (b, d, f) tan/3 = 50. The notations are described in the text. The regions 

outside the inner dash-dotted curves have problematic high-scale tachyons for higher values of M^,. 
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outer curves correspond to = Mp/\/^ and Mqut, respectively. The areas outside 
these contours may be problematic, depending on the cosmological scenario, as discussed 
previously. In panel (a) of Fig. [3l we have chosen e = 1 and tan f3 = 10. The problematic 
region is when |yu| ^ 1000 GeV, reducing to \fi\ ^ 750 GeV for small m^. When e is decreased 
to 10^'^ (panel c) and 10^^ (panel e), the problematic regions extend down to smaller values 
of |yu|, reaching as low as ~ 550 GeV for e = 10~^ and small rriA, essentially independent of 
M,. 

The most immediately noticeable features of the (/i, m^i) planes for tan f3 = 50, shown 
in panels (b), (d) and (f) of Fig. [3l are the greater extent of the stau LSP region when 
fi ^ 1000 GeV, and the the greater extent of the sneutrino LSP region when ^ ^ 1000 GeV. 
In between, the problematic tachyonic regions depend more sensitively on the value of 
than was the case for tan j3 = 10, and also vary more as e is reduced. Once again, it is the 
regions of large /i that are problematic. 

Whereas in the CMSSM case the stau was always lighter than the lightest neutralino 
in the tachyonic region, so that it could be cosmologically acceptable only if the gravitino 
were the LSP, in the NUHM case the tachyonic region also includes parts of the WMAP 
strips where the LSP is the neutralino and it has an acceptable relic density. The crossover 
strip, the stau coannihilation region and parts of the snu coannihilation strip and parts of 
the rapid-anihilation funnel with an acceptable neutralino relic density are all in the non- 
tachyonic parts of the (/i, m^) planes for tan/? = 10. These regions also have acceptable 
b —>■ s'-f for fi > 0. In the case of tan/5 = 50, parts of the stau coannihilation strip and the 
rapid-annihilation funnel are again tachyon-free. We emphasize yet again that the tachyonic 
regions at larger are not necessarily excluded: that would depend on the cosmological 
scenario and whether it avoids the high-scale vacuum in the early Universe. 

In Fig. m we show more examples of {fi, m^) planes, this time with mi/2 = 500 GeV, 
mo = 300 GeV, and Aq = 0. We see that there are regions at larger m^, particularly for 
fi > 0, where the LEP Higgs constraint is respected. As in Fig. [3], the problematic tachyonic 
regions extend to smaller |/i| as rriA decreases, and as e decreases from 1 to 10~^. The 
potential tachyonic problem is restricted essentially to |yu| ^ 1000 GeV. 

Turning to panels (b), (d) and (f) in Fig.|H for tan/5 = 50, we seethat the stau LSP region 
extends to between fi ~ 400 GeV (for small m^) and fi ~ 1000 GeV (for large m^i). The 
problematic tachyon region is now only at larger /i, varying between a range > 2000 GeV 
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FIG. 4: Portions of the NUHM {n,mA) planes for for = 500 GeV, mo = 300 GeV, Aq = 
and (a, c, e) tan/3 = 10, (b, d, f) tan/3 = 50. The notations are described in the text. The regions 
outside the inner dash-dotted curves have problematic high-scale tachyons for higher values of M^. 
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for e = 1 and large tua to ~ 1150 GeV for e = 10~^ and small tua- In all cases, the portion 
of the WMAP strip where the relic neutralino density is controlled by stau coannihilation is 
in the safe region, as well as a portion of the rapid-annihilation funnel where ~ 
Whether the other regions are acceptable would depend on the cosmological scenario. Note 
that the entire regions shown in panels (b), (d) and (f) are favoured by — 2. 

Comparing Figs. [3] and HJ we see (unsurprisingly) that the problematic tachyonic NUHM 
regions grow as mo decreases. As we noted above, it is also possible to consider mg < in the 
NUHM. For example. For mi/2 = 300 GeV and mg/^m^j = -lOO GeV, with = Mqut 
and e = 1, CCB vacua appear at high scale when \^\ > 1000 GeV for rriA ~ 800 GeV 
and > 1500 GeV for lower itla- When e = 10^'^ for the same case, the entire plane is 
problematic. Similarly, when M^, = Mp for any value of e, the entire plane is problematic. 

VI. SUMMARY 

We have discussed in this paper the constraints that are imposed on super symmetric 
models by the presence of the Universe in our familiar EW vacuum. We have argued that 
models with tachyonic spin-zero fields at some high input scale are not necessarily excluded. 
The renormalization-group evolution of tachyonic masses to low scales may change their 
signs, in which case the standard EW vacuum would be a local minimum of the effective 
potential. However, in addition to this vacuum, there may be a much deeper high-scale 
vacuum located along some F- and D-flat direction in the effective potential, with field values 
fixed by some higher-order non-renormalizable interaction. We discuss the circumstances 
under which the existence and location of such a high-scale vacuum can be calculated reliably. 

Such high-scale vacua usually break both color and charge conservation. In general, the 
lifetime for decay of the EW vacuum to this unacceptable lower minimum of the effective 
potential is much longer than the age of the Universe, so future decay into such a vacuum 
is not of immediate concern. A more relevant question is whether the Universe would have 
fallen into such a vacuum during its past history. This depends whether the effective scalar 
masses-squared acquired large positive or negative contributions oc during inflation. If 
these contributions were 0{1) and positive, only the EW vacuum would be populated. On 
the other hand if these contributions were negative (and possibly if they were positive but 
small), the high-scale vacuum would be populated. 
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We have then explored the conditions under which the CMSSM or the NUHM (with 
its two extra parameters) has a calculable high-scale vacuum. If these conditions are not 
satisfied, there is no reason to be tachyophobic. Even if these conditions are satisfied, and 
there is a calculable high-scale vacuum, whether it is catastrophic or not depends on early 
cosmology, and there is still no need to be tachyophobic. 

Unreasoning tachyophobia is never justified: one must examine rationally whether any 
specific tachyonic spin-zero field is dangerous, depending on the evolution of the Universe 
within one's favoured cosmological scenario. 



APPENDIX A: GENERALITIES 



Here we give a concise review of the well-known nature of the fiat directions that will be 
considered. This Appendix is included to serve as a reminder and to set our notation. 

We begin our discussion by ignoring the superpotential and soft terms, which are included 
later. We assume complex scalar fields 0i with a canonical Kahler potential, with vevs 
denoted by Vi. We take Vi to correspond to a supersymmetric vacuum, — 0, in the 
standard notation. Viewed as a vector space, the Vi are the null vectors of the (hermitian) 
mass-squared matrix: 

j:<'^Mi,=j:,i——n=j:,iD'—=o. (ad 

k k,a ^ a * 

Correspondingly, there exists a projection operator into this null space: 

which may be used to construct the modulus field 77, which is the scalar field tangential to 
the null space defined by P: 

^ = E J^2 ^r (A3) 

In addition, there are the other modes Xai that are orthogonal to 77 in the space of the 0j, 
i.e.: 

Xa = Y.'^t<l>-r- Y.^tvi^^. Y.^tdf^S^''. (A4) 
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We now consider the effect of the soft masses on the modulus field 77. When the change of 
variables flA3l) is made, we obtain: 



In the case of a flat direction characterized by a single monomial 

$=nc (A6) 



the D-flatness constraint 



— - • const V i (A7) 



yields: 



Vi\ = Hi ■ const V i. (A8) 



It follows that in this case the modulus soft mass-squared appearing in (1A5I1 is: 

where we see that the weights are nothing but the powers in the monomial. Because we 
specialize to the case of monomial fiat directions flA6p in our study, we use flA9p to determine 
the soft masses of moduli. The cases that we are interested in are those with 

m^{Q) < 0, (AlO) 

where Q is the running scale. In such a case, the modulus runs away from the origin along 
the fiat direction. 

Next we introduce the non-renormalizable superpotential term that lifts the flat direction 
at large field values and stabilizes the modulus against the runaway behavior: 

'*-™-s^,"-E»<.x» + ;^"" (All) 

a 

Here Sa^t are coupling strengths. The fields Xa represent the non-modulus modes (lA4p . 
Below, we make the simplifying assumption that only a single flat direction is 'turned on', 
so that {xa) = 0. 

In Appendix |B] we provide examples that yield the two types of terms appearing in (lAlip . 
We remark that manifest gauge invariance is typically lost when we use the basis r^, Xa, as 
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reflected in the appearance of tlie t]"- coupling. As illustrated by the examples provided in 
the Appendix [Bl the power law dependences on and M^, exhibited in (12.31) are quite 
generic. 

The non-renormalizable superpotential flAlip may well be the result of the exchanges of 
states with mass scale that have been integrated out. The scale M^, could be as low 
as the GUT scale, if appropriate GUT representations are coupled to the MSSM content. 
Alternatively, M^, could be as high as the Planck scale, if it is due to the exchange of 
quantum-gravitational excitations. Another possibility is that ^ Mg 4 to 5 x 10^''' 
GeV, the perturbative heterotic string scale. 

Whatever happens to be, the scalar potential for the modulus r] is given by (12.21) . 
One finds that \u\'^ = |tp + l'^"!^ combines the coupling constants appearing in flAlip . 
Since < 0, the minimum is obtained at: 



w 

In the numerical analysis of flat directions that we perform, we specify the vev v of the 
modulus T] according to the power law that has just been obtained, namely (12.31) . As we 
have commented above, this estimate is consistent with a detailed analysis of the non- 
renormalizable interactions that are allowed in the MSSM. For our purposes, such an order 
of magnitude estimate is sufficient . 



|np(n 



2(ri-2) 



(A12) 



APPENDIX B: EXAMPLES OF FLAT DIRECTION LIFTS 



1. Lift of the HuHd flat direction 

The leading non-renormalizable superpotential term that achieves this is 

Wnr = -^{H^H.f, (Bl) 

where t is a coupling constant. We choose the variant of the fiat direction for which the 
neutral components and get vevs that are equal. In that case the modulus rj and 
orthogonal modes Xa are: 

^ = l^{Hl + H^), xi = ^iH'a-H'd), X2 = H:, X3 = H^- (B2) 
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Using these redefinitions, we find that: 

WNR = ^f + 0{x'), (B3) 

where symmetries forbid XaV^ terms. One such symmetry is ^ H^, which imposes 
symmetry under Xi ^ ~Xi with 7], X2, Xs invariant, which forbids a XiV^ term. The other 
symmetry is —>■ —H^, Hj^ —>■ — which imposes symmetry under rj —1], Xi —Xi 
with X2, Xs invariant. This symmetry forbids X2V^ XsV^- Thus we see that the coefficients 
Sa in ( lAllI) vanish for the present case. The formula ( 1A12I) apphes for the vev, with u ^ t 
and n = 4. Note that fh^ = + {'ni'jj^ + m^^)/2 < must hold for this vev to run away, 
where /i is as usual the (tree-level) MSSM higgsino mixing term. 



2. Lift of the LiL2t'^ flat direction 

We choose the direction where the vevs of 

= ^ = T^ = v, u^ = e = 0. (B4) 

The leading non-renormalizable superpotential term that lifts this flat direction is 

Wnr = -^{HMiLiL2)T', (B5) 

where the parentheses indicate the SU(2) invariants and the coefficient has been selected to 
provide a simple final answer, as becomes evident shortly. As usual, s is a dimensionless 
coupling constant. The modulus r] and orthogonal modes Xa can be parameterized as: 

with X4,5,6 corresponding to the other superfields u^, e, H^. One finds: 

WNR = ^X3v' + Oix')- (B7) 

Thus the stabilization term is of the form of the first term in (lAlip . The rj^ term, with 
coefficient t, is forbidden by matter parity. Symmetry arguments also forbid XaV^ with 
a = 1,2, 4, 5, 6. The formula (lA12p applies for the vev, with u —' s and n = 5. 
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APPENDIX C: METHOD OF SCAN 



In this appendix we summarize the workflow used for our analysis of flat directions. 

1. Starting from unification-scale boundary conditions mo, M1/2, Aq, and electroweak 
scale tan f3, we evolve the RGEs to the EW scale, iterating /i and Bfi until valid 
electroweak symmetry breaking (EWSB) is obtained. We exclude automatically any 
model with tachyonic squarks or sleptons at the low scale. 



2. Next, we loop through all flat directions enumerated in ij] (limiting ourselves to 
monomials flA6|) ) and perform an analysis over the range of running scales Q > 10 
TeV (see the discussion in Section |V]) : 

(a) We check whether < for the weighted sum (lA9p . The powers and masses 

that appear in this sum depend on the flat direction that is chosen; because 
the mass parameters depend on Q, so does the weighted sum m^; 

(b) For values of Q such that < 0, we determine the corresponding vev v us- 
ing (E3D; 



(c) We then compare v to eQ for various values of e: 

• U V < eQ for all Q, our analysis does not show the existence of a CCB 
minimum; 

• Otherwise, the CCB vacuum does exist since the RG-improved analysis tells 
us that there is a potential minimum at large v and the one-loop corrections 
are small over some range of Q for which the vev is nonzero. 
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